Submitted to ApJL 



Self- Absorption of Ionizing Radiation and Extended Narrow-Line 

Emission in High-Redshift QSOs 

S. M. Khairul Alam 1 & Jordi Miralda-Escude 1,2 
ABSTRACT 

We calculate the neutral hydrogen column density of self-absorption in QSOs 
predicted in a model where the QSOs are located in the same halos that contain 
the gas in damped Lya absorption systems. The model is parameterized by the 
probability Pq that any halo has an active QSO. We assume that the QSOs ionize 
the gas, but do not expel or heat it. The derived Hi column densities produce 
negligible Lyman limit absorption, even in the lowest luminosity QSOs, with an 
optical depth of only ~ 10% for luminosity L = 0.01L*, when P = 10 -2 . We also 
compute the Hen Lyman limit self-absorption, which is slightly higher but still 
negligible. The self-absorption can be higher if the gas is highly clumped; only in 
this case the overall emissivity from QSOs could be significantly reduced due to 
absorption by the known damped Lya systems, to affect the predicted intensity 
of the ionizing background or the epoch of Hen reionization. The presence of 
the gas associated with damped absorption systems around QSOs could also 
be detected from the narrow Lya emission line, which should have and angular 
extent of 0.1 to 1" in typical high-redshift QSOs. 

Subject headings: intergalactic medium — large-scale structure of universe - 
quasars: absorption lines 
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1. INTRODUCTION 

The origin of the ionizing background at high redshift has been a long-standing question 
(Bechtold et al. 1987; Miralda-Escude & Ostriker 1990; Songaila, Cowie, & Lilly 1990; Madau 
1991; Haardt & Madau 1996; Giroux & Shapiro 1996). The possible sources include QSO's, 
star-forming galaxies, and the cooling radiation from hot gas in halos. One of the interests in 
measuring the intensity of the ionizing background is that, with the present understanding 
of the Lya forest as arising from the gravitational collapse of structure (see Rauch 1998 for 
a review), the observed mean flux decrement in the Lja forest provides a measurement of 
the parameter where is the baryon density in units of the critical density and T 

is the photoionization rate due to the background (e.g., Rauch et al. 1997, Weinberg et al. 
1997, 1999; McDonald et al. 2000). Once T is known independently, the Lja forest provides 
a measurement of f2& at low redshift which should agree with the values derived from Big 
Bang nucleosynthesis (O'Meara et al. 2001, Pettini & Bowen 2001) and from the Cosmic 
Microwave Background spectrum of temperature fluctuations (Netterfield et al. 2001, Pryke 
et al. 2001, Stompor et al. 2001), if our cosmological ideas are correct. 

The proximity effect (see Scott et al. 2000 and references therein), consisting of the 
reduction of the number of Lja absorption lines in the vicinity of a QSO due to its own 
ionizing radiation, has been used to measure the intensity of the ionizing background, yielding 
a value r ~ 2 x 10~ 12 . This method is subject to several possible systematic errors, including 
uncertainties in the QSO redshift, the effects of gravitational lensing of the QSO (which 
makes it seem more luminous), QSO variability over the photoionization timescale, or the 
clustering of gas around QSOs which may partially balance the ionization effect. 

The simple counting of QSO's as a function of flux in the sky can give us the emis- 
sivity of ionizing radiation, and by taking into account the absorption by the intergalactic 
medium, which is also directly determined in QSO spectra, we can calculate a lower limit 
to the background intensity under the assumption that sources other than QSO's are not 
important (assuming that the escape fraction of ionizing photons from star-forming galaxies 
is negligible). This approach yields T > 10~ 12 , which implies a lower limit VL^h 2 > 0.02 
(Rauch et al. 1997; McDonald et al. 2000, 2001). 

In this paper, we examine if the ionizing radiation from QSOs could be significantly 
absorbed by hydrogen in the halo where the QSO is located, which normally produces the 
damped absorption systems, but is highly photoionized in the presence of the QSO. Although 
absorption by Lyman limit systems is taken into account when computing the intensity of 
the ionizing background, (Haardt & Madau 1996 and references therein), the absorbers have 
always been assumed to be uncorrelated with the QSOs, so any se//-absorption arising in 
the halo of the QSO itself should be in addition to the one computed from the general 
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intergalactic gas. Of course, any such absorption should be directly observable in the QSO 
spectra; however, the self-absorption might be strong only in low-luminosity QSOs, where 
the gas is less highly photoionized, and surveys of Lyman limit systems have generally been 
done on the most luminous QSOs. 

High redshift QSOs are also thought to be the sources that reionize Hen in the inter- 
galactic medium, which could have occurred as late as z = 3 (see Heap et al. 2000 and 
references therein). It is also of interest to know if Hen ions in the QSO halos can produce 
significant self-absorption of the Hen -ionizing radiation from QSOs to delay the epoch at 
which the He n reionization is completed. 

Another interesting consequence of the photoionization of gas in the halo where the QSO 
resides is the extended narrow-line emission that should be produced by the recombinations, 
as recently discussed by Haiman & Rees (2001). We will examine the predicted flux in Lja 
emission. 



2. MODEL 

Our model for the self-absorption in QSO's has two main parts. First, each QSO of 
luminosity L is assumed to be located in a Cold Dark Matter (CDM) halo of mass M, with 
a unique relation M(L), and a probability Pq that any halo will host a QSO. In other words, 
each halo can be either in an active state, in which case the QSO luminosity depends only on 
the halo mass, or in a quiescent state in which case there is no QSO (in reality we expect a 
dispersion in this M(L) relation; we will comment later on its effect on our results). Second, 
a spherical model of the gas density as a function of radius is adopted for each halo of mass 
M, which reproduces the observations of damped Lya systems when the gas is all neutral. 
We then calculate the degree of ionization of the gas in photoionization equilibrium in the 
presence of the QSO flux, assuming that the density has not changed as a result of the 
photoionization heating or hydrodynamic winds from the QSO. 

For the first part of the model, we can determine the relation M(L) once we know the 
distributions of both M and L. We adopt the Press-Schecter formalism (Press & Schecter 
1974; Bond et al. 1991) for the distribution of halo masses, with the CDM model with 
cosmological constant with parameters — 0.7, Q m = 0.3, = 0.04, and h = 0.7. We use 
the parameter 8 C = 1.69 (with top-hat filter) for the threshold overdensity required to form 
a halo, and the relations M = (47r/3) 187T 2 p crit rl ir , and Vj 2 = GM/r vir to relate halo mass, 
virial radius, and circular velocity. 
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For the QSO luminosity function (LF), we use the double power-law model of Pei (1995): 

i{LB;Z) = (L/L B ,y, + (L/L B ,y» • (1) 

where L B is the B-band luminosity. All our results will be presented at z — 3, when 
L B * = 1.2 x 1O 13 L B0 , and 0* = 619.25 GpcT 3 (from Table 1 in Pei 1995, after correcting 
to our adopted cosmological model and Hubble constant). To convert to the luminosity per 
unit frequency at the Lyman limit, we convert to AB magnitude (see Oke 1974) using Mb = 
5.4-2.5 \og{L B /L BQ ) = M AS (4400A) + 0.12 (Schmidt, Schneider, k Gunn 1995), and we use 
a spectral index a = 0.5 from A = 4400 A to A = 1216 A (where L v oc v~ a ), and a = —1.77 
from A = 1216 A to A = 912 A (Zheng et al. 1997). This gives L v * = 1.18 x 10 31 ergs" 1 Hz" 1 
at the Lyman limit. 

Figure 1 shows the relation obtained between the B-band QSO luminosity and the halo 
mass, for four different values of Pq, by simply requiring that the total number of halos with 
mass greater than M times P is equal to the total number of QSOs with luminosity greater 
than L. 

For the density profile of the gas in each halo, we use the same model as McDonald & 
Miralda-Escude (1999): a spherically symmetric gas distribution with an exponential density 
profile, 

p(r) = p exp(-r/r g ) . (2) 

The two parameters of the mass distribution, the radius r g and the density normalization 
p , need to be chosen to match the observed column density distribution of the damped Lya 
systems, when the density p is assumed to be all neutral gas. We use the same parameters as 
McDonald & Miralda-Escude (1999) at z = 3 (see their Fig.3): c g = r g /r vir = 0.04, and p 
determined by a fraction of the baryon mass in the halo in atomic gas form fni = 0.1, both 
assumed to be independent of halo mass. The resulting cumulative number of absorbers 
above a column density Nhi per unit redshift is shown in Figure 2 as the thin line. The 
points in this figure are the observations from Storrie-Lombardi et al. (1996); to obtain these 
points, we divided the total number of absorbers shown in their Figure 4 by their total 
redshift pathlength, Az = 74.7 (which we obtain by adding the redshift pathlength for their 
systems between z = 2 and z — 3, equal to 31.3, and for their systems at z > 3, which is 
43.5, according to their Table 4). 

The calculation of the self-absorption column density is then done simply by calculating 
the neutral fraction as a function of radius for each halo mass M and corresponding QSO 
luminosity L, with the density profile in equation (2), computing the ionizing flux from 
the QSO at each radius, and assuming photoionization equilibrium with the recombination 
coefficient «a = 4x 10~ 13 cm 3 s~ 1 (at an assumed gas temperature T = 10 4 K). 
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3. RESULTS 

The predicted neutral hydrogen column density as a function of the QSO luminosity is 
shown in Figure 3. The Lyman limit optical depth, tll, reaches unity at a column density 
N HI = 1.6 x 10 17 cm -2 , indicated by the thin solid line. In general, t ll is predicted to be quite 
small. Even for L = 0.01L*, t ll ~ 0.1 for P = 1CT 2 . The predicted Hen column density is 
shown in Figure 4. We have assumed a mean spectral index L v oc z/~ L5 between the ioniza- 
tion edges of Hi and Hen , which implies N He ii = 13AN H i (we do not include self-shielding 
effects). Although the Hen Lyman limit opacity is higher than for Hi , it still does not pro- 
duce a very significant absorption. Integrating the quantity L e~ TLL over the QSO luminosity 
function, we find that the overall reduction in the emissivity is a factor (0.85, 0.94, 0.97, 0.99) 
for Hi , and (0.71,0.86,0.93,0.97) for Hen , for P = 10~\ 10~ 2 , 10~ 3 , 10~ 4 , respectively. 
Notice that the reduction in the total ionization rate from the ionizing background is smaller 
than these factors because absorption is lower at frequencies above the Lyman limit. 

It is easy to see that in our model, the Hi column density is proportional to the halo 
mass divided by the QSO luminosity. At fixed r/r g , the gas density is constant, and the flux 
is proportional to L/r 2 , so the neutral fraction goes as r 2 /L, and the column density goes as 
Tg/L. Since r 3 oc r% ir oc M (where r vir is the virial radius of the halo), we have Nhi oc M/L. 
If the halo properties were independent of QSO luminosity, then of course N HI oc Lr 1 . Our 
curves show a slower decrease of the column density with luminosity because of the increasing 
halo mass. The column density also decreases with Pq at fixed L proportionally to the halo 
mass, which is shown in Figure 1. 

Lyman limit absorption at the QSO redshift which decreases with luminosity. A possible 
difficulty is that, in photoionization equilibrium, an equal number of recombinations and 
photoionizations will take place, and about 38% of the recombinations are direct to the 
ground state and produce photons just above the Lyman limit frequency, smoothing the 
discontinuity due to the absorption. As we shall see below, the angular size of the damped 
systems from which these recombination photons would come is probably in the range 0.1 — 1", 
and therefore difficult to resolve from the ground. 

4. DISCUSSION 

We have computed both the neutral hydrogen and helium column densities of self- 
absorption in QSOs in a model where the QSOs are located in the same halos that produce the 
damped Lya absorption systems. The model assumes that each QSO luminosity corresponds 
to a halo mass, and that the gas in damped Lya systems that is ordinarily present in a halo 
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in the absence of a QSO is only photoionized, but not expelled when the QSO is present. 

The parameter that we vary in our model is the probability P that a given halo has 
an active QSO. Large values of P imply long-lived QSOs. For a radiative efficiency e, the 
Salpeter time for the growth of the black hole mass is Mc 2 / (eL E dd) >4x 10 8 e years. The 
typical time for the QSO luminosity function to evolve is ~ 10 9 years; if we require that 
the black hole mass function does not evolve on a shorter timescale, and that e < 0.1, then 
Pq < 0.04. From Figure 3, this limits the self-absorption from the gas in damped Lya 
systems to tll < 0.02 at L*, and t ll < 0.2 at 0.01L*. 

An alternative to the assumption we made of a fixed relation between QSO luminosity 
and halo mass is that, even if QSOs are located in halos, there is little correlation between 
luminosity and halo mass. In this case, we can define a minimum halo mass M min that hosts 
QSOs above a certain luminosity L min , with probability Pq; then, Figures 3 and 4 still give 
the HI and Hell column densities expected at L min , but at higher luminosities the average 
column density would decrease as Nhi oc L -1 , more rapidly than in Figures 3 and 4. 

We note here a possible caveat of our model when P > 10~ 2 : the halos that account 
for most of the observed damped Lya systems have velocity dispersions in the range 40 to 
150 km s -1 , or, at z — 3, total masses of 10 10 to 1O 12 M (e.g., Gardner et al. 1997). From 
Figure 1, this mass range corresponds to rather low QSO luminosities if P ~ 0.01. This 
means that the absorption in QSOs of luminosities > 0.1L* is in this case due to halos with 
masses as large as 10 13 M Q , which could have different physical conditions than the more 
numerous lower-mass halos that are mostly responsible for damped Lya absorbers. 

The Lyman limit self-absorption shown in Figures 3 and 4 could be substantially in- 
creased if the absorbing gas is clumpy. We have assumed that the gas density in the halo is 
smoothly distributed, following the density profile in equation (2). However, the observations 
of multiple metal absorption lines associated with damped Lya systems (e.g., Prochaska & 
Wolfe 1997) show that the gas is actually clumpy. The typical number of multiple absorbers 
that are observed tells us that the covering factor of these clumps is of order unity, but 
the clumping factor depends on the size of the clumps. The clumping factor may be not 
much larger than unity in a scenario like that proposed by Haehnelt, Steinmetz, & Rauch 
(1998), where the clumps are due to halo mergers and the complicated line structure arise 
in a continuous medium of halo gas from velocity caustics and moderate density variations. 
But the clumping factor could be high if the clumps were much smaller than the overall size 
of the damped Lya systems. The absorbing column density increases proportionally to the 
clumping factor. 

Another possible modification of the model we have used here is that the gas observed 



-7- 



in damped Lja systems does not exist in every halo, but only in a certain fraction / of halos. 
In order to preserve the rate of incidence of the observed absorption systems, we then need 
to increase the radius r g in equation (2) by a factor f" 1 ^ 2 , and decrease the central density 
by f 1 ^ 2 to have the same column density distribution in the absence of ionization. With the 
same QSO luminosity, the HI column densities in the presence of the QSO are increased by 
a factor f^ 1 ^ 2 . Thus, making the damped Lya absorbers bigger and less abundant would 
increase the amount of self-absorption, if QSOs are located in the halos that contain the 
damped Lya systems. This could be the case, for example, if QSO activity usually takes 
place in recently merged halos with a lot of fresh gas. 



4.1. Narrow-line emission 

The absorption of the continuum ionizing radiation of a faint QSO by gas with velocity 
dispersions similar to those in damped Lja systems implies that narrow Lja emission lines 
should be observed, containing about a third of the energy that is absorbed in Lyman limit 
photons. This emission line could remain unresolved, since damped Lja systems are small. 
For example, in a halo of mass 10 12 M , r g = 3.15 Kpc (see eq. [2]) at z — 3, corresponding 
to an angular size of 0.41". This narrow emission line would be superposed with the common 
broad absorption lines in luminous QSOs. 

The Lja emission from halo gas around QSOs was recently considered by Haiman & 
Rees (2001). Their model assumes that all the baryons in a CDM halo are in an extended 
gas halo containing a hot and cold phase, with a large fraction of the baryons in the cold 
phase in most halos of interest (see their Fig. 1). In our model, with our choice of parameters 
fui — 0.1 and r g /r vir = 0.04, only 10% of the baryons are in the gaseous halo, but these 
are mostly concentrated within ~ 10% of r V i r . Since the gas density profile in the model of 
Haiman & Rees is roughly isothermal, their prediction for the abundance of damped Lja 
systems in the absence of a central source should not be very different from that of our 
model, which we have shown agrees with the observations. 

Generally, narrow Lja emission may arise from gas over a wide range of radius and 
densities, not only associated with the damped absorption systems but with lower column 
density absorbers as well. The exponential gas density profile in equation (2) for our model is 
intended to approximate the effects of self-shielding when computing the number of damped 
systems in the absence of a central photoionizing source (when the photoionization is due to 
the external background), which should cause the neutral density to drop sharply outside the 
self-shielded region. But the total gas density is likely to drop more slowly. As an example, 
if p g oc r~ 2 , and the gas is mostly ionized, then the recombination rate per unit volume is 
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proportional to r~ 4 , and the Lya surface brightness drops as r~ 3 . This " Lya fuzz" should 
always be present in all halos as long as the gas is ionized, whether or not a central QSO is 
present. In the absence of a central QSO, the Lya emission should have a core at the radius 
where the gas becomes self-shielded against the external background, or in other words, the 
radius at which Lyman limit absorption would be seen against a background source. When 
the QSO turns on, the gas within this core becomes ionized, and the steep power-law surface 
brightness profile can be extended inwards, making it much brighter. However, since Lyman 
limit systems are only ~ 10 times more abundant than damped Lya systems, their typical 
extent in halos should be only ~ 3 times larger in radius than the damped systems, which 
as mentioned earlier would be hardly resolved from the ground. Therefore, the presence of 
a QSO can only increase the emitted Lya in a central region of size ~ 1", except in very 
massive and gas rich halos where the self-shielded region might be larger. 

The Lya surface brightness predicted by our model can be estimated from the H i absorption 
column densities shown in Figure 3. Let us take, for example, the curves for Pq = 10~ 2 in 
Figures 1 and 3, and consider the case of a halo with a = 150kms _1 , which implies a 
mass 1O 12 M at z — 3. Our model associates this halo with a QSO of Lb — 0.03L* ~ 
10 29 5 erg s -1 Hz -1 , which has a Lyman limit absorption tll — 0.04 from Figure 3. The 
frequency over which t ll decreases is Au ~ 10 15 Hz, so the power absorbed in hydrogen 
photoionizations is about 10 43 erg s" 1 , and the power radiated in Lya photons is 10 42 - 5 erg s _1 , 
which at z = 3 produces a Lya flux ~ 10~ 16 erg s -1 cm -2 . If this flux is extended with a profile 
falling as r~ 3 , with a core of ~ 0.5 arc seconds (corresponding to the radius r g of the damped 
absorption system), the surface brightness would be (10~ 17 5 , 10" 19 5 erg s -1 cm -2 arcsec -2 at 
angular separations of (0.5, 3)", respectively. Just like the Hi absorption column density, this 
Lya surface brightness can be increased by the clumping factor of the gas. In comparison, 
from Figure 2 of Haiman & Rees, they conclude that a similar halo (with T vir = 2 x 10 6 K) 
would produce a Lya surface brightness of 2 x 10~ 16 Hz at 3". The large difference with our 
prediction of the Lya surface brightness can be traced to the high clumping factor of the 
cold gas in the model of Haiman & Rees (see their eq. 1), and to their assumption that the 
Lya flux comes mostly from a large region, close to the virial radius of the halo, whereas we 
assume it is much more concentrated to the center. 

Recently, Steidel et al. (2000) have reported the discovery of two "Lya blobs" of emis- 
sion, with surface brightness ~ 10~ 17 ergs -1 cm~ 2 , and angular scale of ~ 10". While the 
surface brightness is similar to the expected value quoted above from a typical damped ab- 
sorber once the gas is ionized, the angular scale is very large. These systems are therefore 
likely produced in exceptionally massive and gas-rich halos. 

To summarize, if extended gas is present in the halos where QSOs are located, with a 
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similar distribution as the average profiles derived from the observed damped Lya systems, 
then we should observe narrow-line Lya emission, extended over a region comparable to the 
size of damped Lya systems, which is typically less than 1". As mentioned at the end of 
§3, there should also be extended emission of the Lyman continuum photons from direct 
recombinations to the ground state. We showed in Figures 3 and 4 that the absorption 
column densities expected are optically thin even for very low QSO luminosities; therefore, 
for fixed halo properties, the Lya fuzz should be easier to observe around fainter QSOs, 
which already produce sufficient flux to ionize all the halo gas, implying a Lya brightness 
independent of QSO luminosity (our Fig. 3 shows that N HI decreases more slowly than L^ 1 
with luminosity, which implies an increasing Lya surface brightness with luminosity, only 
because of the assumption in our model that brighter QSOs are located in more massive 
halos which contain more gas). If the absorption or the associated Lya emission are not 
observed at the predicted level, the conclusion should be that either QSOs are in halos that 
do not contain the average amount of gas that is inferred from the abundance of absorption 
systems, or the QSOs themselves have expelled this gas in winds. 

We would like to acknowledge Patrick McDonald for many stimulating discussions. 
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Fig. 2. — Cumulative number of absorbers per unit redshift as a function of Hi column 
density. Squares are the observed values from Storrie-Lombardi et al. (1996); the thin solid 
line is the prediction of our model. 
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Fig. 3. — Neutral hydrogen column density as a function of QSO B-band luminosity for 
different probabilities Pq. The thin horizontal line indicates a Lyman limit optical depth of 
unity. 
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Fig. 4. — Hen column density as a function of QSO B-band luminosity for different proba- 
bilities Pq. The thin horizontal line indicates a Lyman limit optical depth of unity. 



